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We ask the question, what are the ideal characteristics of a near-held thermophotovoltaic cell? 
Our search leads us to a reformulation of near-held radiative heat transfer in terms of the joint den¬ 
sity of electronic states of the emitter-absorber pair in the thermophotovoltaic system. This form 
reveals that semiconducting materials with narrowband absorption spectra are critical to the energy 
conversion efficiency. This essential feature is unavailable in conventional bulk semiconductor cells 
but can be obtained using low dimensional materials. Our results show that the presence of matched 
van Hove singularities resulting from quantum-conhnement in the emitter and absorber of a ther¬ 
mophotovoltaic cell boosts both the magnitude and spectral selectivity of radiative heat transfer; 
dramatically improving energy conversion efficiency. We provide a model near-held thermophoto¬ 
voltaic system design making use of this idea by employing the van Hove singularities present in 
carbon nanotubes. Shockley-Queisser analysis shows that the predicted heat transfer characteris¬ 
tics of this model device are fundamentally better than existing thermophotovoltaic designs. Our 
work paves the way for the use of quantum dots, quantum wells, two-dimensional semiconductors, 
semiconductor nanowires and carbon nanotubes as future materials for thermophotovoltaic cells. 


I. INTRODUCTION 

The thermophotovoltaic (TPV) method generalizes 
the concept of conventional photovoltaics by replacing 
the source solar spectrum with the thermal emission of an 
engineered selective emitter [T]. When heated, the emitter 
is designed to deliver photonic power over a narrow spec¬ 
tral window which can be efficiently converted into elec¬ 
trical power by a matched photovoltaic cell, assumed to 
be near room temperature. In this way, electrical power 
can be extracted from any heat reservoir [2j [3]. In princi¬ 
ple, the TPV approach avoids the two primary loss mech¬ 
anisms of single-junction solar photovoltaics: the gener¬ 
ation of sub-bandgap photons, and the thermalization of 
excitons with energy much greater than the bandgap[4]. 
For this reason, the theoretical power conversion limita¬ 
tions of TP Vs are near that of an ideal heat engine [5l|6]. 

However, in practice, implementation of the TPV idea 
is difficult. The spectral irradiance of any emitter in the 
far-held is bounded by that of a blackbody with equal 
temperature [7], see Fig.l. For an emitter with tempera¬ 
ture below approximately 1000 K, the fraction of emitted 
power in the range of contemporary low bandgap photo¬ 
voltaics (0.4 — 0.75 eV) is negligible [8], making practical 
photovoltaic conversion problematic. To overcome this 
limitation, the emitter of current TPV systems must op¬ 
erate at a temperature between roughly 1200 and 1500 
K, creating severe design requirements. Specihcally, the 
emitter must be thermally robust in this temperature 
range, suppress emission below the bandgap of the pho¬ 
tovoltaic cell, which constitutes the majority of the solar 
spectrum, and have emissivity near the blackbody limit 
above the bandgap pfTO]. Even without considering tem¬ 
perature induced degradation of optical properties and 
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structure, current theoretical proposals still require emit¬ 
ter temperatures surpassing 1500 K to reach 40% conver¬ 
sion efficiency under Shockley-Queisser analysis pThlS] . 

An intriguing prospect to partially alleviate this high 
emitter temperature requirement is to bring the matched 
emitter and absorber into the near-held; where the pres¬ 
ence of tunneling photons allows heat transfer to sur¬ 
pass the far-held blackbody limit (near-held thermopho- 
tovoltaics, NFTPVUMS]). Making use of matched 
near-held resonances pTl [TS] . for example surface plasmon 
polaritons[T9l [20] as in Fig.l, radiative heat transfer can 
be made spectrally sharp and signihcantly greater than 
the far-held limit (Super-Planckian emission [21]). 

Nevertheless, basic design issues persist even in theo¬ 
retical consideration. Analysis of near-held radiative heat 
transfer has conventionally focused on the maximization 
of total transfered power, leading naturally to designs us¬ 
ing coupled modes of symmetric structures. However, to 
extract usable electrical power, one half of the near-held 
resonant architecture must function as a photovoltaic 
cell. And when realistic bulk semiconductors are included 
as the low temperature absorber, radiative heat transfer 
is greatly reduced and spectrally broadened [22] [23]. For 
example, considering a 100 nm gap, and identical hypo¬ 
thetical Drude metals with surface plasmon polariton res¬ 
onances matched to the bandgap of gallium antimonide 
photovoltaic cell [24], replacing the symmetric metal ab¬ 
sorber with the photovoltaic cell reduces radiative heat 
transfer by a factor of nearly 10, Fig.l. Furthermore, al¬ 
though the temperature diherence between the emitter 
and the photovoltaic cell can be comparatively reduced, 
bringing the photovoltaic cell into the near-held of the 
emitter will, in practice, increases its operating tempera¬ 
ture and thus reduces its conversion efficiency [25]. Resul- 
tantly, much of the performance gain offered by moving 
to the near field thermophotovoltaic systems is lost. 

In this article, we show that the narrowband optical 
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FIG. 1. (A) Comparison of radiative heat transfer between 
the near-held and far-held calculated using 0- The near-held 
system consists of two half-spaces with optical response of a 
hypothetical Drude metal, one at 750 K, the other at 0 K, with 
0.08 eV loss parameter, 0.58 eV plasma frequency, and an ad¬ 
ditional static imaginary permittivity of 1, separated by a 100 
nm gap. This arrangement is shown schematically as an inset 
in panel (B). The vertical gray lines mark an approximate 
range of current single-junction low bandgap photovoltaics. 
The enhancement of radiative heat transfer through evanes¬ 
cent modes relaxes the necessity of maintaining the emitter at 
extremely high temperature. (B) Comparison of the enhance¬ 
ment of radiative heat transfer beyond the blackbody limit for 
half-spaces separated by a 100 nm gap, shown schematically 
as an inset. The enhancement is temperature independent. 
The dashed curve is found by assuming identical hypotheti¬ 
cal Drude metals with 0.21 eV loss parameter, 0.83 eV plasma 
frequency, and an additional static imaginary permittivity of 
2.5. The solid curve results when one of the metal half-spaces 
is replaced by a gallium antimonide photovoltaic cell. The 
inclusion of the bulk photovoltaic cell greatly reduces the en¬ 
hancement effects provided by operating in the near-field. 


absorption arising from quantum confinement leads to 
spectrally selective and efficient near-field heat transfer, 
see Fig.2. This is a fundamental advantage since ab¬ 
sorbed thermal power is the source of exciton genera¬ 
tion. Simultaneously, we also uncover the root cause 
of the poor performance of traditional bulk semicon¬ 
ductor photovoltaic cells in the near-field by recasting 
near-field radiative heat transfer in terms of the joint 
density of electronic states. Appealingly, this results 
requires only simple arguments from Shockley-Queisser 
analysis and the essential Kramers-Kronig optical re¬ 
sponse constraints. We also provide the basis for an ex¬ 
perimentally realizable near-field thermophotovaltic sys¬ 
tem making use of carbon nanotubes, and calculate its 
enhanced performance metrics. Our predictions can be 
verified by near-field thermal emission spectroscopy, pro¬ 
viding a road map for the use of quantum dots, quan¬ 
tum wells, two-dimensional semiconductors, semiconduc¬ 
tor nanowires and carbon nanotubes as the future build¬ 
ing blocks of thermophotovoltaic devices. 


Selective Emitter: 



Electronic Absorber: 




FIG. 2. Schematic illustration of the central result of this 
article. Here, ED OS stands for the joint density of electronic 
states. We will show that by matching peaks in the joint 
density of electronic states of the emitter to peaks in the joint 
density of electronic states of the photovoltaic cell the spectral 
selectivity of near-field radiative heat transfer and efficiency 
of near-field thermophotovoltaics can be greatly improved. 


II. ENERGY CONSIDERATIONS FOR 
MAXIMAL RADIATIVE HEAT TRANSFER 


From the pioneering works of Rvtov[26]. and Polder 
and van Hove [27]. the spectral irradiance between two 
planar half spaces separated by a vacuum gap is given by 
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Where c is the speed of light in vacuum, 0 {Ti^ujp) the 
mean energy expectation value of the canonical harmonic 
oscillator at temperature and angular frequency ujp^ d 
the vacuum gap separation normalized by the magnitude 
of the free space wave vector at uOp^ Ti the temperature 
of the ith half space, kp the parallel component of the 
wave vector, again normalized by the magnitude of the 
free space wave vector at ujp^ the j-polarized reflec¬ 
tion coefficient of the ith half space, k\_ the normalized 
perpendicular component of the wave vector in the vac¬ 
uum gap, and Im[...] the imaginary part of the enclosed 
function. 

Intuitively, radiative heat transfer between planar half 
spaces is maximized when the structures are symmetric, 
and the vacuum gap is made vanishingly small, d —> 0. 
In this limit the integrand of the spectral irradiance. 
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HiT,,T 2 ,LUo,kp) from Eq.Q, is dominated by the high 
momentum p-polarized evanescent contribution, /Cp ^ 1 , 
and is well approximated by the bound form 


H {Ti,T2,ujo,kp) 


aj^,\e{T^,ajo))-e{T2,aJo)\ _ 

47r2c2 ^1 + (e'(wo) /e"(wo))^^ 


( 2 ) 

with e (ujo) denoting the relative permittivity of both me¬ 
dia, and the ' and " superscripts marking the real and 
imaginary parts of the function. To achieve the greatest 
possible radiative heat transfer between identical media 
at given wave vector and frequency, the so called upper 
bound of radiative heat transfer per channel [28]. the ratio 
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must be made as small as possible. (This condition is 
antithetical to the plasmonic field enhancement figure of 
merit where 7 {uq) is maximized.) 

Applying the Kramers-Kronig relation [29j to the real 
part of the permittivity, the 7 { 00 o) factor can be ex¬ 
pressed entirely in terms of e (ujo) as 
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For energy harvesting with semiconductor photovoltaic 
cells, Shockley-Queisser efficiency analysis can be used 
to provide constraints to the material parameters which 
minimize j(ujo)' (1) No radiative heat transfer should 
occur below the bandgap frequency cjg because it cannot 
be converted into useful electrical power. Mathematically 
this criterion is stated as j (cUo) cxo for all cUo < 
equivalent to e {ujo) 0 for all uJq < ojg. The dispersive 
part of 7 {ujg) then provides an explicitly positive con¬ 
tribution, and its minimization requires e {uJo) 0 for 
all uJo 7 ^ ujg. This restriction on the global absorption 
characteristics of the media, e (cjo), can be immediately 
translated into two statements concerning its polariza¬ 
tion, e{uJo)‘ First, e{ujg) > 0. Second, the polarization 
of the media acts as a store for photonic power. From 
this store the initially transferred power may either be 
converted into the internal degrees of freedom of the ab¬ 
sorber or return to the emitter, limiting the total trans¬ 
ferred power. The larger the polarization of the medium 
is the smaller radiative heat transfer will be. 

( 2 ) radiative heat transfer should be made as great as 
possible at the bandgap frequency, e (ujg) 00 . This 
second requirement for achieving 7 {ujg) ^ 0 is provided 
by the fluctuation dissipation theorem. To maximize ra¬ 
diative heat transfer at we require e {ujg) ^ oc to 
generate the largest possible thermal currents, the first 
term of Eq. The balance of these two conditions leads 
unequivocally to the conclusion that to minimize 7 
under the above Shockley-Queisser constraints, e {uJo) 


must be sharply resonant about ujg. Strictly, simultane¬ 
ous minimization of 7 { 00 g) and complete suppression of 
heat transfer for all sub-bandgap frequencies is only pos¬ 
sible if e {uJo) is mathematically equivalent to the Dirac 
delta distribution, e {(Jq) ^ i^g)- 

The implications of this result to near-held energy har¬ 
vesting become immediate by recalling the intimate con¬ 
nection between optical dissipation, and the joint density 
of electronic states. 
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with ao denoting the polarization vector of the elec¬ 
tromagnetic excitation, pji the matrix element of the 
momentum operator, cOq the angular frequency, £ji (k) 
the constant energy surface between hlled and excited 
states such that £ji (k) = Tiujo^ dS^ the inhnitesimal 
constant energy surface element in d dimensions, and e 
and rrie the mass and charge of an electron. Applying 
Eq.([^ to model free particle semiconductor systems in 
one, two and three spatial dimension we recover the well 
known energy scaling behaviors e (uOo) oc ^JT^uJo — Eg {d= 
3), C (d = 2), and 1/ — Eg{d = 1 ), where (7 is a 
constant, and Eg the energy of the bandgap, match¬ 
ing the behavior of the joint density of electronic states. 
It is important to note that suppression/enhancement 
of electronic states at a given frequency is the key to 
decreasing/increasing optical absorption. However, the 
Kramers-Kronig relation introduces dependence on the 
global frequency characteristics through the polarization 
of the medium. Both these factors have to be taken into 
account in order to minimize the factor governing spec¬ 
trally selective near-field heat transfer, 

As previously mentioned, the principal method for 
overcoming the Shockley-Queisser limit in TPV devices 
is to spectrally tailor the photonic output of the emit¬ 
ter to frequencies just above the bandgap of the photo¬ 
voltaic cell. But, in the near-field this frequency range 
also corresponds to the poorest absorption characteristics 
for a typical bulk semiconductor. Near the bandgap the 
— Eg scaling of the joint density of electronic states 
gives rise to both low absorption, e and large po¬ 

larization, which can be seen by the positive contribution 
of principal value part in the Kramers-Kronig relation¬ 
ship (cjo) = 1 + f /o^ duo ^ • These two fac¬ 

tors lead to large 7 {uJq) factors detrimental to spectrally 
selective heat transfer when uOo is near uOg^ Eig.l. In short, 
the usual properties of a three dimensional semiconduc¬ 
tor are in complete disagreement with those required for 
maximal radiative heat transfer as given by Eqn. Q. 

In stark contrast, significant suppression of the ^{uOq) 
factor can be obtained if the absorption, e (cjo), of the 
semiconductor cell is strongly peaked at the operating 
frequency. According to Eq. such optical behavior 
is achieved when the joint density of electronic states 
becomes similar to the ideal Dirac distribution at the 
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bandgap. Suppression of the electronic states below 
the operating frequency is necessary to decrease trans¬ 
fer of inefficient sub-bandgap photons, as dictated by 
Shockley-Queisser analysis. The suppression of states 
above the bandgap is necessary to decrease the polar¬ 
ization of the medium at the bandgap, as required by 
the Kramers-Kronig relations. These features can essen¬ 
tially be obtained through van Hove singularities, which 
occur in a free particle model in low spatial dimensions. 
This evolution is shown in Fig.3, where the 7 {ujo) fac¬ 
tors for free particle semiconductor systems in one, two 
and three dimensions have been calculated. This struc¬ 
ture in the density of electronic states and narrowband 
absorption behavior is characteristic of semiconductors 
with quantum-confined spatial dimensions such as quan¬ 
tum dots, nanotubes, and two-dimensional materials. We 
conclude based on the above analysis that systems pos¬ 
sessing these properties are ideally suited for achieving 
spectrally selective radiative heat transfer and can func¬ 
tion as ideal near-field thermophotovoltaic cells. 



FIG. 3. : Base ten logarithm of the longitudinal 7 {uo) 
e {(jUo) /e (<^o) factor for one, two, and three dimensional 
model free particle semiconductors. The inset shows a 
schematic illustration of the evolution of the joint density 
of electronic states for zero, one, and two quantum-confined 
spatial dimensions characteristic of a semiconductor based on 
the work of Kiimmel et al[30]. By introducing dimensional 
constraints the dispersion of the joint density of electronic 
states can be altered to produce spectrally selectivity super- 
Planckian radiative heat transfer near the bandgap of a photo¬ 
voltaic cell. This result illustrates both the primary drawback 
of employing traditional bulk semiconductor photovoltaics for 
near-field photonic energy conversion, and the usefulness of 
switching to photovoltaics with quantum-confined dimensions 
for this application. 


III. MODEL THERMOPHOTOVOLTAIC 
SYSTEM 

A wide body of previous work exists in the areas of 
photovoltaic cells that utilize materials such as quantum 


dots, graphene and carbon nanotubes [3Tti33] . Our anal¬ 
ysis in the previous section shows that any of these low 
dimensional material systems are excellent candidates for 
future near-field TPV cells. Again we emphasize that the 
efficiency boost arises from spectrally selective nature of 
the absorption in these systems which is the hallmark of 
quantum confinement. 

As a framework of such a device, we consider here 
an idealized NFTPV system consisting of an emitter- 
absorber pair utilizing carbon nanotubes [34] (CNTs) and 
calculate the radiative heat transfer characteristics which 
can be observed in experiment. Our CNT model device is 
analogous to the one-dimensional free particle semicon¬ 
ductor example shown in Fig.3. The confinement of elec¬ 
trons in CNTs leads to strong van Hove singularities [35]. 
spectrally selective absorption, and semiconductor be¬ 
havior for specific chiral vectors. Note that nearly identi¬ 
cal performance would occur if the CNTs were exchanged 
for semiconductor nanowires [36]. We have focused on 
this CNT design, despite current challenges in fabrica¬ 
tion and electron-hole pair collection [37], primarily be¬ 
cause the thermal robustness [22138| and low bandgap of 
CNTs seems ideally suited to TPV applications. 

Two versions of the CNT system have been consid¬ 
ered. The first is designed for a 1300 K emitter and 300 
K absorber; the second, for an 800 K emitter and 300 
K absorber. In both cases, the absorber is composed 
of CNTs with chiral vector (19,0) embedded vertically 
in a matrix of zirconium dioxide on a tungsten backing. 
This choice of chiral vector allows the CNTs to operate 
as semiconductors, with a bandgap of 0.58 eV. This be¬ 
havior occurs as the Dirac point of the band-structure is 
avoided [ 39 ] due to the angular quantization. The exact 
current voltage characteristics of the CNT photovoltaic 
cells have not been included, and instead ideal p-n junc¬ 
tion behavior has been assumed. This is in accordance 
with the Shockley-Quiesser analysis for ideal performance 
limitations. (To create fully functioning NFTPV cells the 
CNT on tungsten absorbers would be attached to exter¬ 
nal electronics and act as the photovoltaic elements [4QF 
02) 

Both emitters are composed of free standing CNTs 
with the same chiral vector as the CNTs considered for 
the absorber, and are set on a zirconium dioxide back¬ 
ing. The absorber and emitter are separated from each 
other by a 100 nm vacuum gap. In the high tempera¬ 
ture case the fill fractions and thickness of the emitter 
and absorber CNT layers are 65% - 50 nm, and 65% - 20 
nm. For the low temperature case these parameters are 
altered to 30% - 50 nm, and 65% - lOOnm in the same 
order. A schematic of this setup is included as an inset 
in Fig.5. 

The effective perpendicular permittivities of the two 
considered designs are shown in Fig.4. The absorption 
spikes, e can be tuned by varying the chiral vector. 

The first principles calculation of these optical properties 
is based on the Kubo formalism described by Falkovsky 
and Varlamov [ 43 ]. with additional loss included via the 
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relaxation time approximation with an estimate relax¬ 
ation time of one picosecond [44l [45]. The effective par¬ 
allel permittivities, which play a secondary role in deter¬ 
mining the transferred power, are calculated in identical 
manner. The main steps of this calculation are outlined 
in the appendix. Again, a variety of near-held emitters 
using other resonances, such as surface plasmon polari- 
tons or phonon polaritons, or absorbers, such as those 
mentioned above, could also be used with the CNT ab¬ 
sorber to obtain very similar results. 
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FIG. 4. : Effective relative permittivity perpendicular to the 
half spaces for the CNT emitter absorber pair designed to 
operate at 1300 K, depicted in Fig.5. For vacuum gaps on 
the order of hundreds of nanometers radiative heat transfer 
is dominated by tunneling photons with p-polarization. In 
this case, the overall radiative heat transfer properties are 
primarily determined by the perpendicular permittivity. As 
predicted by Eqn.(2), the radiative heat transfer is maximized 
when the 7 (c^o) factor is minimized for both the emitter and 
absorber. 

In hgure 5, we plot the spectral irradiance for this CNT 
based emitter-absorber system calculated using Rytov’s 
theory. As can be seen by comparing the perpendicular 
permittivity. Fig.4, with the spectral irradiance. Fig.5, 
the spectral position of peak radiative heat transfer cor¬ 
responds strongly with the van Hove singularities and 
the minimization of j(uJo)‘ In fact, although the ana¬ 
lytical theory of is valid only in the very near-held, 
the spectral position of peak photonic transfer varies only 
4% as the vacuum gap is increased up to roughly half the 
wavelength of operation, 1 jum. 

The spectrally tailored nature of emission and ab¬ 
sorption can lead to signihcant enhancement for energy 
conversion. Applying Shockley-Queisser analysis [4] as 
an upper-bound performance estimate, the CNT system 
with the 1300 K emitter (Fig. 5) could supply up to 15.00 
W/cm^ of electrical power with 53.0% power conversion 
efficiency. Likewise, the 800 K emitter system could sup¬ 
ply up to 1.07 W/cm^ of electrical power with 15.6% 
power conversion efficiency. For comparison, an opti¬ 
mized conventional bulk gallium antimonide[24] NFTPV 
design has also been considered. This system consists 
of a 175 nm thick Drude metal, with 0.74 eV plasma 
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FIG. 5. : radiative heat transfer for two metamaterial emit¬ 
ter absorber pairs utilizing semiconducting GNTs optimized 
for emitter temperatures of 1300 K and 800 K. In both in¬ 
stances the absorber is assumed to be held at a temperature 
of 300 K. Parameters for the two systems are given in the 
main text. A schematic of the system is inset over the plot. 
The greatly enhanced, spectrally thin, radiative heat trans¬ 
fer peak produced slightly above the energy bandgap of the 
matched semi-conducting GNT absorber seen here is ideal for 
near-field photonic energy conversion applications. 

frequency, 0.21 eV loss parameter, on a tungsten back¬ 
ing with a vacuum gap separation of 100 nm. Under 
identical Shockley-Queisser analysis at a 1300 K emitter 
temperature this conventional system could produce 1.62 
W/cm^ of electrical power with 18.8% power conversion 
efficiency, whereas with an emitter temperature of 800 K 
it could create only 0.01 W/cm^ of electrical power with 
1.6% power conversion efficiency. More concretely, 73.1% 
of photonic thermal power transferred from the emitter 
in the CNT NFTPV device lies above the bandgap, 20.7 
W/cm^, for the 1300 K model system and 25.0%, 1.7 
W/cm^, in the 800 K model. In comparison, the opti¬ 
mized plasmonic system described above transfers 28.1% 
of its thermal power above the bandgap of the gallium an- 
timonide photovoltaic cell, 2.4 W/cm^, when the emitter 
is given a temperature of 1300 K and 3.2%, 0.02 W/cm^, 
when the emitter is assumed to have a temperature of 
800 K. 


IV. SUMMARY 


We comment briefly on previous far-held TPV designs 
which have utilized Gallium Antimonide [Ml EH as the 
photovoltaic cell for converting the thermal radiation into 
electric power [48] . Our analysis shows that the optimum 
cell design for near-held TPV is fundamentally diherent 
and the presence of van Hove Singularities in the material 
comprising the cell is critical for the spectrally selective 
nature of the transferred energy. This in turn should lead 
to fundamental improvements in energy conversion effi¬ 
ciency as shown by the Shockley Queisser analysis. The 
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implementation of near-field TPV designs is more chal¬ 
lenging than far-held TPV designs. However, a funda¬ 
mental promise of near-held TPV, as mentioned before, 
is enhanced heat transfer for lower operating tempera¬ 
tures. This arises because the black body limit which 
fundamentally constraints far-held TPV does not apply 
to near-held TPV where heat exchange takes place due 
to tunneling of evanescent waves. 

In summary, we have shown that the ideal joint den¬ 
sity of electronic states for near-held photonic energy 
capture is mathematically equivalent to the Dirac delta 
distribution. This result immediately reveals why bulk 
semiconductors are clearly ill suited to near-held pho¬ 
tonic energy capture. Thus switching to photovoltaic 
cells with van Hove singularities, seen in any semicon¬ 
ductor with a quantum-conhned dimension, ohers a clear 
path for improving the efficiency of future NFTPV de¬ 
vices. For experimental verihcation, we have provided 
model designs for two such low dimensional systems con¬ 
sisting of matched metamaterial CNT emitter absorber 
pairs and have shown that their Shockley-Queisser limit 
performance metrics are well beyond those achievable 
with current NFTPV device designs using bulk semi¬ 
conductor photovoltaics. Our main aim is to emphasize 
that the absorption characteristics of such low dimen¬ 
sional materials are ideal for near-held TPV and is an 
interesting avenue to explore in the development of ther¬ 
mophot ovoltaic cells. Our work will also motivate fu¬ 
ture studies to adapt quantum well, quantum dot, two- 
dimensional semiconductors and CNT cells specihcally 
for thermophotovoltaics. 
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VI. APPENDIX: OUTLINE FOR THE 
CALCULATION OF THE OPTICAL 
PROPERTIES OF CNTS 


As mentioned in the main text, the calculation of 
the optical properties of the (19,0) chiral vector CNTs 
used in our sample NFTPV designs has followed the 
Kubo formalism method described by Falkovsky and 
Varlamov [43]. The central equation for the dielectric per¬ 
mittivity tensor in this approach is 

n \- , Pab{k,uj) 

^ab \k^UJ) ^sur T 7 : (h) 

hceo UJ 


where 


p.a«<.k) = / 

Tr [Va (p) g K,p+) Vb (p) g K,p-)] 


(7) 


Here, Fab denotes the polarization matrix, with the sub¬ 
scripts a and b indicating directions, Q the Matsubara 
Green function, uj the frequency of the exciting photon, 
ujn the discrete frequencies of the Matsubara sum, Tr the 
trace operation, the permittivity of free space, Csur the 
relative permittivity of the surrounding medium, Va the 
velocity operator matrix defined by the Heisenberg oper¬ 
ator evolution equation 




Xa = 


[Xa,H] ^ dH 
ih dpa ’ 


( 8 ) 


p the momentum integration variable, and the + and — 
superscripts the values p + | and p — | with k standing 
for the momentum of the photon excitation. Moreover, 
the Green’s function is explicitly defined as 

g = {iwn-H{p))-\ (9) 


and its the sum over the discrete frequencies of the Mat¬ 
subara sum, in terms of its matrix elements indexed by 
the subscripts i and j, by 


E 


Qi Qj ) 


fo jp )) - fo (p^)) 

UJ F ip - {Sj (p^) - £i (p-))' 
( 10 ) 

In these definitions H is the Hamiltonian of the system, 77 
the loss parameter following the relaxation time approx¬ 
imation, and fo {£j) the Fermi distribution at the energy 
level of the j eigen state. 

For the carbon nanotube system treated in the text, 
only TT-orbital interactions have been considered due to 
the energy range of interest [49]. With these approxima¬ 
tions the Hamiltonian is 


. . / -ipaVS inj • pa , iTvj P^\ 

H{p,j)=^[e 2 n g Vs j ^ (11) 

where a is the lattice constant for graphene, 2.46 A, j the 
angular number, p the momentum of the electron along 
the tube, n a constant related to the wrapping of the 
tube|34], and 7 is the orbital overlap energy for graphene, 
3.1 eV. The loss parameter 7 = h/r is estimated with a 
relaxation of one picosecond [44[ |45] for r. The energy of 
the nanotube is similarly defined as 


£j = ±7 


\ 


1 + 4 cos 


( \fZpa 



+ 4 



( 12 ) 

To extend the theory presented by Falkovsky and Var¬ 
lamov from graphene to the single layer nanotubes, the 
Hamilitonian associated with graphene is first enlarged 
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to take into account the energy level splitting resulting 
from the additional confinement in the nanotube. This 
step is accomplished by the substitution 




0 

0 

0 

'OH' 


0 

0 

to 

0 

H* 0 


0 iJ|2 0 0 



. 0 0 0 . 


for the axial direction, and 



- 0 

0 

ves2 

0" 


"0 

0 

0 

0 - 

Ce+ = 

0 

^982 

0 

0 

0 

0 

0 

0 

Vs- = 

0 

0 

0 

0 

0 

0 

"^981 

0 


0 

0 

0 

0. 


.0 

^9 El 

0 

0 


(15) 

for increasing and decreasing angular momentum along 
the angular direction respectively, see [34]. Note that for 
the angular velocity operators the derivative is taken with 
respect to the natural discrete momenta. These operators 
are then converted to the eigenbasis of the Hamilionian 
by the linear transformation v with 


with H as above, and the £i and £2 subscripts denoting 
the lowest level energy bands for the particular chiral 
vector chosen following the discussion provided by Wong 
and Akinwadej39]. The velocity operators for axial, z, 
and angular direction, 0, are then: 


U = 


Hei 

\Hei\ 

0 

0 

1 


0 0 

H£2 He2 

\H£2\ \H£2\ 

1 1 

0 0 


Hei -I 
\Hei\ 

0 

0 

1 


(16) 


An identical procedure is undertaken for the effective 
mass operators, Mij = dp dp- ^ calculation 

of the static conductivities used in the relaxation time 
approximation. With these substitutions the procedure 
detailed by Falkovsky and Varlamov [43]. equations ^ 
and can be followed directly to achieve Fig.4. 


0 0 0 vzsi ' 

0 0 VzS2 0 

0 0 0 

0 0 0 . 


(14) 
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